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Abstract 

Baryons made of three heavy quarks become weakly coupled, when all the quarks are sufficiently 
heavy such that the typical momentum transfer is much larger than Aqcd- We use variational 
method to estimate masses of the lowest-lying bcc, ccc, bbb and bbc states by assuming they are 
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CN ' Coulomb bound states. Our predictions for these states are systematically lower than those made 
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I. INTRODUCTION 



The field of heavy quark spectroscopy is experiencing a rapid renaissance, mainly pro- 
pelled by the emergencies of several unusual charmonium resonances, of which X(3872), 
F(4260) are the highlights Accompanied with these unexpected discoveries, progress 
has also been made steadily in the more traditional sector of charmonium spectroscopy, 
exemplified by the recent sightings of several long-awaited particles such as rj c (2S), h c , and 
particularly the doubly- charmed baryons such as H+, Precise knowledge of their prop- 

erties will help to refine our present understanding of heavy quark dynamics Q]. 

After the tentative establishment of the doubly charmed baryons |3|] , one may naturally 
expects to fill the baryon family with the last missing member, i.e., baryons composed 
entirely of heavy quarks, denoted the QQQ states in short. Being a baryonic analogue of 
heavy quarkonium, the triply- heavy baryons are of considerable theoretical interest, since 
they are free of light quark contamination and may serve as a clean probe to the interplay 
between perturbative and nonperturbative QCD. 

One of the basic properties of these heaviest baryons in Nature is their masses, which will 
be the primary concern of this paper. In contrast with the spectra of the singly-heavy and 
doubly-heavy baryons, to which a vast number of literature based on either phenomenological 
approaches or lattice QCD simulations are dedicated, only sparse attention has been paid 
to the spectroscopy of triply-heavy baryons, perhaps mainly due to the lack of experimental 
incentive. 

The interest toward these baryons can be traced back to Bjorken, who first carried out 
a comprehensive studies on their properties two decades ago, particularly focusing on the 
discovery potential of the triply-charmed baryon state Reconstructing a QQQ candidate 
is a rather challenging job experimentally, since it is difficult to separate all the decay 
products emerging from the cascade decay chain QQQ — ► QQq — > Qqq from the copious 
hadronic background. Nevertheless, according to Bjorken, some semileptonic decay channels 
such as — > + 3/i + + 3z/ M , may offer a clean signature for a ccc event. 

Needless to say, the discovery potential of triply-heavy baryons also crucially depends 
on the production environment. Baranov and Slad have shown that the production cross 
sections for triply-charmed baryons at e + e~ collider are too tiny to be practically relevant j^. 
Gomshi-Nobary and Sepahvand have recently calculated the fragmentation functions of c 
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and b evolving into various triply-heavy baryons, and estimated that the corresponding 

n 

fragmentation probabilities vary in the range 10 7 ~ 10 4 6]. They consequently estimated 
two largest cross sections, which are associated with producing Qf, cc and fi ccc , to be about 
2 and 0.3 nb in the forthcoming Large Hadron Collider (LHC) experiment with cuts of 
Pt > 10 GeV and \y\ < 1. For an integrated luminosity of 300 fb _1 (about one year of 
running at the LHC design luminosity C = 10 34 cm 2 s _1 ), the amount of flbcc and Q ccc yield 
can reach about 6 x 10 8 and 1 x 10 8 . It seems rather promising to establish these two states 
in such a large data sample. 

Stimulated by the discovery possibility of triply-heavy baryons in near future, it is no 
longer of only academic interest to study their properties like mass spectra. Unfortunately, 
no predictions to the masses of triply heavy baryons from lattice QCD simulations have 
emerged yet (only the static three-quark potential has been measured ^0]), and one has 
to resort to other theoretical means at this moment. 

Heavy quarkonium spectroscopy is traditionally the arena of phenomenological potential 
models, which in general incorporate a long-range confinement interaction jslll0j,lu| (see also 
Godfrey's contribution in Q]). Nevertheless, recent advances in nonrelativistic effective field 
theories of QCD, particularly the effective theory dubbed potential NRQCD (pNRQCD), 
has started to put h eavy quarkonium spectroscopy on a model-independent ground [12j (for 
a recent review, see |l3j). A novel aspect of this effective field theory is that the interquark 
potential arises as the matching coefficients. In this language, different quarkonium states are 
categorized with respect to the relative magnitude between the typical momentum scale, mv, 
and the nonperturbative QCD scale, Aqcd- In the case of mv ^> Aqcd, the corresponding 
state is said to be weakly coupled, and the dynamics is largely dictated by the short-distance 
potential which can be calculated order by order in a s ; in the other situation like mv ~ 
Aqcd, the state is said to be strongly coupled since the potential is no longer calculable 
in perturbation theory, instead must be determined by nonperturbative methods such as 
lattice QCD. It is in the latter situation that a pNRQCD framework intimately resembles 
the phenomenological Cornell model. Evidences are accumulating to hint that T, B c may 
well be identified with the weakly-coupled system, whereas J/ ip lies in the borderline between 
the weak and strong couping regime, and the first few excited bottomonium states (far from 
the open flavor threshold) belong to the strongly-coupled system j"!^ . 

In parallel with the formulation of pNRQCD for the quarkonium, Brambilla, Vairo and 
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Rosch have recently laid down an analogous framework for triply-heavy baryons [14|. The 
effective Lagrangian has been written down for both weakly-coupled and strongly-coupled 
QQQ states, with some of the matching coefficients supplied. Among various possible appli- 
cations of Ref. Q], exploring the mass spectra of QQQ states is the most straightforward to 
think of. This is the very goal of the present work. To make things more tractable, we will 
confine ourselves in this work to the weakly-coupled states only. The Qttt, if exists, would 
be an ideal prototype for such a state. However, to be phenomeno logically relevant, we have 
to stick to baryons made exclusively of bottom or charm. As in quarkonium, most probably 
only the ground states are amenable to a weak-coupling assignment. To be objective, due 
to weaker interquark color strength in a baryon than in a meson, and not so heavy charm 
and bottom masses, one cannot exclude the possibility that even the ground states might 
be strongly coupled. 

Despite this disclaimer, we will proceed by assuming that the QQQ ground states are 
indeed the weakly-coupled system. In this work, we attempt to estimate the leading order 
contribution to the binding energy, therefore for this purpose, only the tree-level static inter- 
quark potential, i.e., Coulomb potential, needs to be considered. Since rigorously solving a 
three-body Coulomb bound state problem is beyond our current ability, we have to resort 
to some sort of approximation method. Stimulated by success of the variational method in 
coping with few-body atomic system, we will invoke this simple but efficient approximation 
scheme to address our baryonic problem. For baryons containing simultaneously b and c, 
we will take advantage of the mass hierarchy nib 3> m c to guide our variational analysis, 
just in analogy with that in the simple 3-body atomic system such as helium atom and the 
ionized hydrogen molecule, the physical picture becomes much more tractable by exploiting 
the fact Tfi 7v 3> m e . 

The rest of the paper is distributed as follows. In Section [TT1 we present a brief introduc- 
tion to the most relevant features of the triply- heavy baryons in the weak-coupling regime. 
In Section ITTTj which is the main body of this work, we perform a detailed variational analysis 
to the binding energy of various QQQ ground states. Three different classes of triply-heavy 
baryons, bcc, ccc (bbb) and bbc states are treated separately, with the hierarchy m;, ^> m c 
utilized as a guidance for choosing proper trial states. Considerable amount of effort has 
been devoted to the bbc state, which is the most interesting case, but also most difficult to 
analyze. We have employed three different approaches to study this state and explored the 
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implication of each approach in depth. Particularly the relevance of the compact diquark 
picture is discussed. In Section IIV1 we present our predictions to the masses of all the 
lowest-lying triply heavy baryons, and compare our results with other work. We summarize 
and present an outlook in Section |3 



II. WEAKLY COUPLED QQQ STATES 

In this section we recapitulate the major aspects of weakly-coupled triply heavy baryons 
which are most relevant to this work. For more comprehensive discussion from the perspec- 
tive of pNRQCD, we refer the interested readers to Ref. 



13, 



14] 



To efficiently investigate the low energy properties of a tripled-heavy baryon, such as 
binding energy, it is convenient to work with a low energy effective theory that focuses 
on the most relevant degrees of freedom. In a weakly-coupled QQQ state, the relevant 
low energy degrees of freedom are nonrelativistic heavy quarks and (ultrasoft) gluons with 
energy and momentum of order mv 2 , just like in a weakly-coupled quarkonium. All the high 
energy degrees of freedom, which can only appear in virtual states, have been integrated out 
explicitly. One particularly important high energy mode is the (soft) gluons with momentum 
of order mv, whose effects are encoded in the low energy theory as the interquark potentials. 
Since we have mv ^> Aq C d m a weakly coupled state, the potentials can be determined in 
perturbation theory by matching procedure. There are infinite number of potentials, which 
are organized in expansions of 1/m. The most important potential is the O(l/m ) (static) 
potential. 

The explicit form of potentials depends on the overall color configuration of quarks. Three 
heavy quarks can be in either color singlet, octet or decuplet state. The color-singlet state 
represents the most important case, since it constitutes the leading Fock component of a 
baryon. The singlet static potential is well known, 

vj 0) = ( ^— + — + ^— ) + 0(*l) , (1) 

3 \Jxi-x 2 | |x 2 -x 3 | |x 3 -xi|y 

where the color interaction between any pair of quarks is attractive. In general, in the color 
octet and decuplet configurations, some or all pairs of quarks will repel each other. 

Thus far, the low energy effective theory is completely depicted by a set of uncoupled 
Schrodinger equations governing the motion of heavy quarks in different color configura- 
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tions. The situation becomes more intriguing when ultrasoft gluons are included. Since the 
typical wavelength of ultrasoft gluons is much longer than the typical interquark distance, 
the gluon fields can be multipole expanded. Very much like the electromagnetic multi- 
pole transition in atoms, ultrasoft gluons can also induce chromo-electromagnetic multipole 
transition from one heavy quark color configuration to a different one. In particular, a 
chromo-electric dipole transition can occur between a color-singlet QQQ configuration to an 
octet one. The interaction between bound heavy quarks and vacuum gluonic fluctuations 
through this chrome- El operator, will generate the leading nonperturbative correction to 
mass of a triply- heavy baryon, as a manifestation of Lamb shift in QCD . The magnitude 
of this nonperturbative correction depends on the relative size between mv 2 and Aqcd- It 
is quite difficult to estimate this effect accurately, and we will not consider it further. 

We end this section by commenting briefly on the solidity of the weak-coupling assignment 
to the lowest-lying triply heavy baryons that are of practical interest. As was admitted in 
Introduction, since the interquark color strength in a baryon is only a half of the quark- 
antiquark color strength in a quarkonium, the typical dimension of a triply heavy baryon, 
say, Qbbb, is expected to be considerably fatter than that of T. One may worry that the 
wave function of the former could permeate deeply into the confinement region. Fortunately, 
as will be shown quantitatively in the forthcoming section, the interquark attraction is 
effectively enhanced due to the influence of the third quark, so the actual situation turns 
out to be considerably better than this pessimistic anticipation. 

III. VARIATIONAL ESTIMATE OF BINDING ENERGY 

In this section, we attempt to estimate the binding energy of various QQQ ground states. 
Our starting point is the color-singlet hamiltonian: 

where the ellipsis stands for the higher-dimensional potentials suppressed by powers of 1/m. 
Because the purpose of this work is to calculate the leading 0(a 2 s ) contribution, we will 
restrict to the lowest order static potential only. 

1 The analogous effect in quarkonium system, originally considered by Voloshin [15j, has been extensively 
explored by many authors. 
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To describe a bound state, we need first separate the relative motions of quarks from 
the center-of-mass motion in (J2J. There are infinite ways to perform this separation. A 
simple way is to replace the old coordinates by the center-of-mass coordinate plus two new 
coordinates defined as the positions of the quark 1, 2 relative to the quark 3: 



m,-x 



ri = xi - x 3 , 

r 2 = x 2 - x 3 . (3) 

In terms of these new coordinates, the hamiltonian (J2J) can be separated into 

H s = H™ + h s , (4) 

where Hg M = -V|/ (2 X] rnj) is the center-of-mass part, and the part governing the relative 
motion reads 

, s = _ZL_^L_Z^_^fI + I + M (5) 

2m 13 2m 23 m 3 3 Vn r 2 r 12 J 

where r 12 = |r x — r 2 |, = (1/mj + l/m,) -1 is the reduced mass between quark i and j. In 
such a coordinate system, the quark 3, sitting at the origin, is artificially singled out from 
two other quarks. 

Our task then becomes solving the bound state problem defined in (jSJ). In the following, 
we will use the variational method to estimate the corresponding binding energy of each 
type of QQQ ground states. 



A. bcc 



We start from the simplest case, the baryon made of one heavier bottom quark of mass 
M and two lighter charm quarks of mass m (throughout this work, we will use the notation 
M = mb and m = m c ). 

It is convenient to choose a coordinate system as specified in (J3J), with b sitting at the 
origin. This coordinate system is sketched in Fig. [T] Subsequently, substituting mi = m 2 = 
m and m$ = M into (jSj), the hamiltonian describing the internal motion of a singlet bbc 
state reads 

2m red 2m red M 3 \r 1 r 2 ) 3 r 12 ' 
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FIG. 1: Sketch of the coordinate system used for the bcc state. 

where m re( j = (1/m + 1/M) -1 is the reduced mass between c and b. Note in this choice of 
coordinates, the motion of b is embodied in the reduced mass and the operator Vi • V2/M. 

Let us first consider an ideal bcc state with M/m — > 00. In this situation, the b quark 
just acts as a static color source, with two c quarks revolving around. This picture is very 
similar to that of the two-electron atoms such as H~ , He and Li + , where the nucleus is 
practically fixed in space, and two K-shell electrons orbit about it. Estimating the energy 
of the two-electron atoms is considered as a classical application of the variational method, 
which has been discussed virtually in every quantum mechanics textbook (e.g., see jl6j|). 

Closely following the textbook treatment of helium, we may approximate the bcc ground 
state to be the one in which each of the c moves in the Is orbital of an effective Coulomb 
potential, somewhat stronger than — 2a s /3r. This is so because the attractive color interac- 
tion felt by each c due to b is strengthened due to the attraction exerted by another c. This 
is in opposite situation to He, where nuclear charge felt by each electron is partly screened 
due to the repulsion exerted by another electron. 

In a physical bcc state, the hierarchy M ^> m is much less perfect than that in a helium. 
Nevertheless, the above ansatz about the form of the ground state wave function still seems 
plausible. What we need is to take the motion of b into account. For notational simplicity, 
we will take the "baryonic" unit m re d = 2a s /3 = 1, in which all the length and energy scales 
are measured in the unit of the Bohr radius (m re d 2a s /3)~ 1 and Bohr energy m re d (2a s /3) 2 . 
We choose the spatial part of trial wave functions as 

*(ri,r 2 ) = /(r 1 )/(r a ), (7) 
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where 

\3/2 

f(r) = ~^e- Xr (8) 

'IT 



is the normalized Is Coulomb wave function. Here A is a variational parameter, which 
characterizes the effective color charge of b perceived by each of the c. Obviously, when the 
attractive interaction between two charm is turned off, A would simply be 1. For the reason 
discussed earlier, we expect A > 1 in our case, so that each c can be thought of moving on 
a squeezed Is orbital. This is opposite to what is expected in a helium. 

In the He ground state, two K-shell electrons must form a spin singlet to obey Fermi 
statistics, since its spatial wave function is symmetric under the interchange of two electrons. 
Due to the extra color degree of freedom carried by quarks, two c quarks in the bcc ground 
state must instead be a spin triplet. When combined with b, the lowest-lying bcc baryon can 
be either J p = | + or J = | + , which are degenerate up to 0{m?a^/M) corrections due to 
the hyperfine splitting. 

We now attempt to find the expression for the ground state energy, E. Taking the 
expectation value of hs with the trial wave function in (JZJ), after some effort one obtains 

E = -X 2 + 2X(X-1) + J, (9) 

where 

J = -~, d 3 n d 3 r 2 e - = -£A, (10) 



7T 2 J J r i2 
measures the average potential energy stored between two charm quarks. 

Note that the double integral involving the Vi • V2 term vanishes, because of spherical 
symmetry possessed by the Is wave functions. Thus, the effect of kinetic energy of b is fully 
taken into account by the reduced mass m re d- 

The minimum of energy can be found by requiring dE/dX = 0, which leads to 

A = — , (11) 
16 ' v ; 

indeed compatible with our expectation. The corresponding ground state energy is 

E = -(§) 2 ^-Q 2m -d« s 2 , (12) 

where the Bohr energy has been inserted in the last term, to recover the actual dimension 
of energy. 
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B. ccc 



The triply charmed baryon states no longer have an atomic counterpart. On the other 
hand, the ccc ground state is highly constrained by symmetry. To have lowest energy, it nec- 
essarily possesses a totally symmetric spatial wave function. After the totally antisymmetric 
color wave function is included, Fermi statistics then demands that it must have J p = | + . 

We again work with a coordinate system defined in Q, with one c, artificially denoted 
charm 3, fixed at the origin. The hamiltonian depicting the relative motion of three identical 
c can be obtained by making the replacement M — > m in (JSj). We then have the reduced 
mass m re( j = m/2. To condense the notation, we will work with the "baryonic" unit, in 
which the corresponding hamiltonian becomes 

We are attempting to seek a proper form for the trial wave function for ccc ground state. 
One simplest choice is motivated from that adopted for a bcc state. Let us temporarily 
imagine the charm 3 can be distinguished from the rest of two, then (JJ|) constitutes a 
reasonable representation for such a state. Now coming back to a physical ccc state, to 
account for the indistinguishablity of c, we should fully symmetrize 0. With the spin part 
of wave function suppressed, the trial wave function then reads 

U fr . x _ /(ri)/(r 2 ) + /(ri)/(r 12 ) + /(r 12 )/(r 2 ) 

*(ri, r 2 ) . = , (,14) 

/3(1 + 2T) 



where / is the same as given in (jSji . and contains a variational parameter A. One can check 
this wave function is symmetrical under the interchange of any two charm quarks. The 
is normalized by incorporating the overlap integral T, 

T = ^ II e- A(2ri+r2+ri2) = gjj • (15) 

Taking the expectation value of hg in the trial state (|14jh after some straightforward 
manipulation, we end up with the expression 

-A 2 + 2A(A - 1) + J - A 2 T + 2(A - 1)V - 4Q + T - Q 

1+2T ' ^ ' 

where J has been given in (JTUJ), and 

V = — d 3 ri d 3 r 2 = —A, 

7T 2 J J T\ 81 
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\6 A(2ri+r2+ri2) 

Q = * d 3 ri d 3 r 2 e 



it 2 J J r 2 

A 6 ff ,o ,o e~ A ( 2ri+r2+ri2 ) 16, 

= — // d i rid i r 2 = — A . 

Wi n 2 27 

.F = [(d'nd^e-^^Vle-^ = -^A 2 , 

7T Z jj 243 

= ^ ffd i r 1 d i r 2 e- x ^ + ^V l -V 2 e- x ^ + ^ = ^-X 2 . (17) 
7T^ jj 243 

Note the exchange integrals V, Q T and arise from the symmetrization effect, which are 
absent in the expression for the energy of the bcc baryon, (jHJ). In particular, Q, the double 
integral involving Vi ■ V2, no longer vanishes this time. 

Substituting the results of these integrals into (fTtjf). we obtain 

2595 531 

952 595 K J 

The optimum can be found by variational principle, 

_ 4325 ^ ^ 527 (19) 



2 832 



and 



3 741125 „ , . 

E = - 17Q7376 — > -0.925 m red a 2 , (20) 

where the normal unit is recovered in the last entity. 

It is interesting to compare the results we have got for the bcc and ccc ground states. 
First lowering the b mass in a bcc state down to m, we get a fictitious ccc state with one 
c distinguishable from the other two. Comparing (J2(J|) and (|12|). we immediately find the 
symmetrization effect tends to lower the energy. Moreover, by comparing (|19|) and (flip, we 
find the symmetrization effect also tends to compress the bound state size more. 

For actual bcc and ccc states, we find E^ ccc > En bcc (note m re d in two cases are different), 
which implies that charm quarks in f2& cc are more tightly bound than in Q ccc . This is con- 
sistent with the general expectation that a bound state with constitutes of vastly disparate 
masses is more stable than that with equal-mass constitutes, say, a hydrogen atom is more 
stable than a positronium. 
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FIG. 2: Sketch of the coordinate system adopted for the bbc state. 



C. bbc 

We finally turn to baryons made of two heavier b quarks and a lighter c quark. This type 
of baryon is more complicated than the preceding two, because the effective potential felt 
by c is no longer spherically-symmetric, but merely axially- symmetric. 

To make the symmetry between two b quarks manifest, we may adopt a more appropriate 
coordinate system other than @. Letting mi = 7712 = M, 1713 = m, we define the following 
new coordinates: 

M(xx + x 2 ) + mx 3 
~~ 2M + m ' 

R = xi - x 2 , 

Xl +x 2 

r = x 3 . (21) 

Note now the coordinate origin coincides with the middle point between two b quarks. The 
geometry of these new coordinates can be clearly visualized in Fig. EJ Substituting (|21|) 
into the original hamiltonian (J2J), we find that the part of hamiltonian responsible for the 
internal motion is 

ftg = __ZL_?a.l__ZL_*!./I + r), (22) 
2M red 3 R 2m red 3 \n r 2 J 



where M red = M/2, m red = (1/m + 1/2M) 1 are the reduced masses, and r\ 



Ir - 

I 1 2 1' 



r2 = |r+ f I 



A nuisance may deserve some caution before we move on further. Two strong coupling 
constants in (|22j) have been tacitly assumed to be evaluated at the same renormalization 
scale /1. This procedure seems incompatible with our intuition that the first a s should be 
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affiliated with a scale ~ 1/(R), and the second one with a different scale ~ l/( r )j where 
(R) and (r) represent the typical values of R and r^, respectively. When M and m are 
widely separated, (R) <C (r) is expected, and this recipe will miss the contributions of large 
logarithm ln((r)/(i?)), no matter which value of /i is chosen. The symptom encountered 
in this equal-a s ansatz is a typical shortcoming of lowest-order perturbative calculation in 
a multi-scale problem, and one in principle can ameliorate its prediction by appealing to 
renormalization group equation to resum the large logarithms of the form a™ ln n ~ 1 ((r) / (R)). 
Fortunately, for a physical bbc baryon, M is only three times larger than m, (R) and (r) 
likely don't differ much, hence we don't need worry much about this nuisance. 

In the following, we will treat the bbc ground state with three different approaches: point- 
like diquark approximation, Born-Oppenheimer approximation, and variational method. 



1. Point-like Diquark Approximation 

In a fictitious world where M is many orders of magnitude heavier than m, the physical 
picture simplifies enormously. The influence of c to the motion of very heavy b can be safely 
neglected. Consequently, two b quarks in the bbc ground state form a Is spin-triplet state. 
The very compact Bohr radius of b, cannot be resolved by c which is orbiting from far away. 
Thus from the perspective of c, the bb cluster is just like a point particle. To distribute 
itself in the lowest energy, the c is revolving around this point particle in the corresponding 
Is orbital. Being in 3 color state, this compact diquark may be identified with a heavy 
antiquark. In this sense, the bbc ground state is analogous to the heavy quarkonium B c . 

In passing, it is worth mentioning that the doubly heavy baryons, such as bbq states, fit 
into this compact diquark picture to a better extent than the bbc state, because the average 
distance between q and diquark in the former, ~ 1/Aqcd, is considerably larger than the 
average distance between c and diquark in the latter, ~ 1/ (ma s ). Properties of doubly-h eavy 
baryons was first studied within a compact diquark picture long ago in HQET language |17| . 
Some refinement to this picture, which invokes the nonrelativistic EFT of QCD to describe 



the internal excitation of the diquark, has recently come out 



14. 



18] 



The form of ()22|) is particularly convenient to accommodate the compact diquark pic- 
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ture 2 . Because (R) <C (r) in this case, one may approximate r\ and r2 by r, the color 
potential felt by c then becomes — 4a s /3r, as if it is due to a heavy antiquark sitting at the 
origin. Eq. (|22j) then collapses into two separate hamiltonians, one governing the internal 
motion of the diquark, the other governing the motion of c in a central Coulomb potential. 
The energy of the bbc ground state is then simply the sum 



For a physical bbc state, the mass hierarchy between b and c is far from ideal, so the 
usefulness of this oversimplified approximation is doubtful. 

2. Born-Oppenheimer Approximation 

We now seek an alternative method that explicitly incorporates the effect of finite diquark 
size. First observe that an ideal bbc state bears some similarities with the simplest molecule, 
the H2 ion, in the sense that both are three-body bound states held together by Coulomb 
force, and both contain two heavy particles and one much lighter particle. Motivated by 
this similarity, one may wonder whether some well-known method developed to analyze 
can be transplanted here. 

A standard tactics to cope with diatomic molecules, such as the ion, is Born- 
Oppenheimer approximation (adiabatic approximation). This method was originally moti- 
vated by the strong separation of time scales between electronic and nuclear motion, which is 
mainly a consequence of the hierarchy m e <mjv. The recipe of this method is that, to solve 
molecular problem, one first determines the electronic eigenstates at fixed nuclear positions, 
then takes the corresponding electronic energy as an effective potential, in conjunction with 
the internuclear Coulomb potential to describe the nuclear motion. 

There is a caveat, however. Despite the aforementioned similarities, one should realize 
there is one fundamental difference between H2 and the bbc state, that is, the internuclear 
Coulomb interaction is repulsive, whereas the Coulomb interaction between b is attractive. 
This difference in turn results in drastically distinct properties of H% and an ideal bbc 

2 For pedagogical purpose, in the following two a s in 12L'I) will be simply taken equal even in the limit 




(23) 



M/m — > 00. 
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state. As a result, success of adiabatic approximation to the former does not automatically 
guarantee that it can be taken for granted for the latter. 

To better orientate ourselves, it is instructive to recall first how an adiabatic picture arises 
from the ground state Q]. A snapshot of this simplest molecule is that, two nuclei slowly 
vibrate about some equilibrium positions with small amplitude, whereas the electron flies 
around much more swiftly. The vibrational nuclear motion is a consequence of the balance 
between internuclear Coulomb repulsion and an effective attractive interaction induced by 
the electron. A crucial fact is that the typical period of nuclear motion is much (~ \JrriN /m e ) 
longer than that of electronic motion. It is thus a good approximation to regard nuclei as 
frozen when considering the electronic motion, consequently the electron will distribute itself 
in the ground state of this static nuclear potential. Moreover, the electron can be regarded 
as responding instantaneously to the change of nuclear arrangement, therefore it follows the 
nuclear motion adiabatically, which implies that it can always remain in the corresponding 
ground state for each nuclear configuration. 

In contrast, an ideal bbc state bears a completely different structure. As we have known, 
this state is characterized by a compact diquark picture. The pull exerted by c again induces 
an effective attractive interaction between two b quarks. However, when superimposed on 
the attractive Coulomb interaction, it helps, though with a rather minor impact, to push 
two b closer. The only agent to prevent a complete collapsing is the kinetic energy of b. It 
is interesting to compare the overwhelmingly dominant role enjoyed by the kinetic energy 
of b with the insignificant role played by the nuclear kinetic energy in H^- 

Based on the point-like diquark picture, one can show that b is confined in a region 
about m/M smaller than c, the typical velocities of b and c are about equal (~ a s ), and 
the typical kinetic energy of b is about Mj m larger than that of c. Obviously, notions such 
as "fast c" and "slow 6" are simply misnomers. Moreover, uncertainty principle tells that 
the typical orbiting period of b is much (~ M/m) shorter than that of c. As a result, c can 
hardly follow the fuzzy pace of b, let alone to readjust itself instantaneously to the ground 
state for a particular configuration of b. In sharp contrast with H2, the bbc state exhibits a 
completely anti- adiabatic nature. 

The above negative argument seems to persuade us to give up adiabatic approximation 
in analyzing an ideal bbc state, since the orthodox picture on which this method is based 
is badly violated. Ironically, this method practically does yield correct result for this state. 
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The reason can be traced as follows. We have argued that it is difficult for c to react quickly 
to the rapid change of configurations of b. It simply gets confused. However, the really 
important point is, what c can see is only a smeared bb cluster which is well localized in a 
small region, it doesn't care about the details going on inside this cluster. What c can do 
is to distribute itself in the ground state of the Coulomb potential due to a remote 3 color 
source. This is of course nothing but the point-like diquark picture. Born-Oppenheimer 
method takes the energy eigenvalue of c in static configurations of b as effective potential for 
b. For an ideal bbc state, only the value of this effective potential at very small separation 
of b is relevant, which is just the Is energy of c in the Coulomb potential of an antiquark. 
Following the Born-Oppenheimer procedure, the motion of b is described by a new potential, 
which is the original Coulomb potential scaled up by this tiny constant. One then readily 
reproduces the correct answer, (J23j) . for the ground state energy. 

It is now clear that Born-Oppenheimer approximation practically works for an ideal bbc 
state because of very compact diquark size. But we certainly are more interested in the 
physical bbc state. Since m and M are not widely separated in this case, there is no more 
strong separation of time scales, this approximation thus is not expected to yield accurate 
result in the first place. Nevertheless, since this method takes the finite diquark size effect 
into consideration, which is relevant for a physical bbc state, we will take a practical attitude, 
applying it to this state to watch what results will come out. 

Let us now concretely analyze the bbc state following Born-Oppenheimer method. To 
start, we first approximate the full wave function \l/ as 



where (p represents the charm ground state for a static configuration of 6, and $ stands for 
the amplitude to find b in this configuration when c is in the state ip. 

In the Born-Oppenheimer ansatz, we need first determine the lowest eigenstate <p by 
solving the Schrodinger equation 



Since the positions of b explicitly enter the potential, the energy eigenvalue e depends on R. 

This is exactly the same problem as one encounters in H2, to determine the electronic 
ground state at fixed nuclear positions, therefore we can follow the standard treatment [19?] . 



*(R,r) « $(R)^(R,r), 



(24) 




(25) 
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FIG. 3: A and effective potential as functions of R, determined by the variational calculus (solid 
line). All the numbers are given in the lighter "baryonic" unit. In the lower half plot, the dashed 
curve represents the function given in (|3U|) . which is hardly distinguishable from the actual one. 

Solving (J2*5|) rigorously is unfeasible, one commonly appeals to variational method. A rea- 
sonable form taken for the trial wave function is a linear combination of Is charm states 
centered on each of b quarks. A variational parameter A is included in the Is trial state to 
characterize the effective color charge of b perceived by c. Taking the indistinguishableness 
of b into account, the trial wave function of c takes the form 

fin) ± /(r 2 ) 



p(R,r) 



(26) 



I2(1±S(X,R)) 

where / is given in (JHJ). To keep our notation simple, we have chosen to work with the 
lighter "baryonic" unit m re d = 2a s /3 = 1. The overlap integral is incorporated to make <p 
normalized, 



S{X,R) 



A ;; 



7T 



d 3 re -X(ri+T2) 



X 2 R 2 



1 + XR + 



-XR 



(27) 



The wave function (p must be either symmetric or antisymmetric upon interchange of 
two b (R — > — R), so that the corresponding bb pair, if in relative s-wave, must be either 
a spin triplet or singlet in line with Fermi statistics. When combined with c, the former 
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configuration corresponds to a bbc state with J p = | or ~ , and the latter corresponds to 
a state with J P = \ . As is well known in H2, the antisymmetric configuration has higher 
energy level than the symmetric one. Since we are only interested in the bbc ground state, 
we will discard the state with odd parity. 

We thus choose the symmetric one in f)26j) . Multiplying both sides of (|2*3j) by the corre- 
sponding tp*, integrating over r, we find that the charm energy reads 

A 2 X(X-1)-C(X,R) + (X-2)S(X,R) 

£{R) = ~y i+s(x,n) • (28) 

The classical interaction integral C and exchange integral S are given by 

\3 , p -2\rt 1 

C(X,R) = - d 3 r- = - l-(l + Ai?) e - 2Ai? , 

7r J r 2 R 1 J 

\3 r p -A(n+r 2 ) 

£(X, R) = — d 3 r = A (1 + XR) e~ XR . (29) 

ix J r 2 

To locate the minimum of (|28|) at a given R, we resort to the condition de/dX\n = 0. The 
analytical expression for the optimum, if can be worked out, would be very cumbersome, so 
we are content with providing numerical solutions only. The optimized A and e as functions 
of R are shown in Fig. El In digression, we would like to mention that a trick adopted by some 
texts (for example, (19J), which aims to facilitate finding the optimum, is mathematically 
inconsistent, therefore we have refrained from using it. 

Fig. El illustrates some expected features of charm ground state in a static configuration 
of b. At R = 0, the bb diquark shrinks to a point, the color charges double, so we have 
A = 2 and e = — 2 2 /2 = —2. This is exactly what we would expect by replacing a point-like 
diquark with an antiquark. As R gets large, c will be essentially localized with one of the 
b, forming a Is state, and the influence of the other b becomes negligible. To put in a 
quantitative way, at large R, the effective charge A ~ 1, and the energy of c is the Is energy 
plus the potential energy between c and the other b, that is, e m — 1/2 — 1/R. 

In the Born-Oppenheimer ansatz, the charm energy plays the role of effective potential 
for b. To expedite our analysis, it is convenient to have an analytic formula that mimics the 
actual e(R), which is known only numerically. We find the following parameterization, 

£& JR) = -0.5 L5 - -.- , (30) 

v ' 1 + 0.586 R 1A21 ' v ; 

represents a good fit to the actual one, with error less than one percent provided that R is 
not too large. As already pointed out, due to the compact diquark nature of a bbc state, 
only the knowledge in small R range affects the bound state property. 
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The remaining task is to determine with the effective potential taken as input. In 
Born-Oppenheimer approximation, the motion of b is simply governed by the following 
Schrodinger equation 3 : 



VI 1 



R 



r> , k < -0.5 - 1.5 
R 



1 + 0.586 (kR) 



1.421 



-1 



$(R) = £$(R), (31) 



where k = m re d/M re d. For convenience, we have switched to the heavier "baryonic" unit 
M re d = 20:^/3 = 1. Note k plays the role of scale conversion factor. 

This equation can be solved numerically once k is specified. Consequently, the energy of 
the baryon ground state, E, can be identified with the eigenvalue of the corresponding Is 
state. The dependence of E on n in a wide range is shown in Fig. HJ As is expected, at small 
k, the energy predicted from this approach does coincide with the one from the point-like 
diquark approximation. Technically, this can be understood by examining (|31j) in the limit 
K — ► 0, in which the effective potential reduces to a constant —Ik. As discussed before, the 
underlying reason should be attributed to the fact that for small k, only the value of the 
effective potential near R = is relevant. 

In short, the lesson we have learned is that, even though Born-Oppenheimer approxima- 
tion is not theoretically justified for an ideal bbc state, this procedure still leads to correct 
results because of the compact diquark nature in this state. 

As k increases, the average diquark size becomes comparable with the typical distance 
between c and b. In this situation, neither point-like diquark approximation nor Born- 
Oppenheimer approximation is expected to make reliable prediction. Nevertheless, since 
the latter approach explicitly incorporates the effect of finite diquark size, we may expect it 
is closer to the truth than the former. As one can discern in Fig. the prediction of E from 
the latter approach becomes incrementally higher than that from the former as k increases. 
This is compatible with our expectation. The larger k is, the more relevant the contribution 
of the effective potential at large separation of b becomes. Since e monotonically increases 
with R (see Fig. |3J), thus Born-Oppenheimer approximation predicts higher E. 



3 To arrive at this formula, one has dropped two additional terms containing V_r ip (see for detailed 
derivation). This procedure is partly justified by the fact V Rip <C V# $ as expected from the compact 
diquark picture. Although a rigorous mathematical proof is absent, we invoke the fact that this procedure 
makes correct prediction for an ideal bbc state as an evidence for its validity. 
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3. One-Step Variational Estimate 



We have shown that both the point-like diquark approximation and Born-Oppenheimer 
approximation render correct predictions for an ideal bbc state. However, there is no a priori 
reason to expect them to work satisfactorily for a physical bbc state, where m and M are 
not so widely separated. A useful indicator is the ratio of the average diquark dimension to 
the typical distance between c and b, which is roughly 

(R) 2Mj,y , . 

±-L ~ K ss — 4^- ~ 0.6. 32 

Because of poor separation between (R) and (r), a more general approach is called for to 
analyze the physical bbc state. 

In the following we will employ the third method, dubbed one-step variational estimate. It 
takes basically the same variational ansatz as used for the bcc and ccc system. However, due 
to more complex nature of the bbc system, two variational parameters have been introduced. 
The term one-step implies that the ground state energy as well as the full wave function are 
determined in a single step, in contrast with Bohr-Oppenheimer procedure, in which one 
determines the wave functions of c and b in two successive steps. On general ground, one 
expects this method is more accurate than the other two, inasmuch as it is based entirely on 
the variational principle and no other approximation is invoked. As long as the trial wave 
function is reasonably chosen, we expect it will render reliable prediction even when « is not 
small. 

For notational convenience, we adopt the heavier "baryonic" unit M re< j = 2a s /3 = 1 here. 
The hamiltonian (}2*2*|) then simplifies to 

V 2 B 1 1 (VI k k\ 

hs = - „ - - hr + - + - , 33 

2 R k \ 2 ri r 2 J 

where the scale conversion factor is included in the c sector. 

We first need to guess a proper form for the trial wave function 1 J r . It is natural to follow 
the ansatz of (|24)l. to express ^ in a quasi-separable form <&(R) <y?(R, r), where $ represents 
the b wave function, and ip denotes the c wave function, which may be taken the same as 
(|2b|). This form of trial wave function clearly embodies the point-like diquark picture in the 
k — > limit. Since ip has incorporated the effects of finite diquark size, this choice of trial 
wave function seems reasonable also for large k. We take the trial wave function for the bbc 



20 



ground state explicitly to be 



1 



£3/2 (ft . A) 3/2 



e -8R ( e -K\ Tl + e -n\r 2 



71 



(34) 



'2(1 + 5) 

with the spin wave function suppressed. A and 5 are variational parameters. Note \I/ is 
symmetric under the reflection R — ► — R, as it should be for the ground state. The wave 
function is normalized by incorporating the overlap integral 

16 5 3 (25 2 + 55nX + 4/t 2 A 2 ) 



IX 



Jd 3 Re- 25R S(K\,R) 



(35) 



(25 + kA) 5 
where S is given in (|27|). 

The physical implication of A is the same as in Born-Oppenheimer ansatz, which describes 
the effective charge of b perceived by c, except here it is taken as a constant instead of a 
function of R. This simplification seems plausible, at least for small k. As noticed before, 
the typical time scale characterizing the change of configurations of b is in general shorter 
than that of c, consequently c only sees smeared trajectories of b. When considering the 
impact of b on c, it is reasonable to average its effects over different configurations of b. This 
averaging procedure will lead to a constant value of A. 

The new parameter, 5, is introduced simultaneously to characterize the impact of c on 
the geometry of the diquark. It would simply equal 1 in the limit k — > 0, when the influence 
of c becomes completely negligible. 

Taking the expectation value of h s , (J3~3~j) . in the trial state 
forward manipulation, we obtain 

k\ 2 k 2 \ 2 



after some straight- 



E=- S - 
2 



+ 



8{8 -!){! + X) 

K 2 \ 



A (A- 1) -C+ (A-2)£] + — (\ + S-48y 



/(I + s). 



(36) 



The parameters in (}3l)j) are given by 

5 



kC = — fd 3 Re- 2SR C(K\,R) 

7T J 



K~E — — — l<l*Rr- 2SH 



E(k\R) 



5k\ (5 2 + 3Sk\ + k 2 X 2 ) 

(5 + kA) 3 
1Q5 3 k\ (5 + 2kX) 



8X 



a 3 

- id 3 Re 



(25 + kXY 

2SR S(kX, R) _ 4 5 3 (45 2 + 85kX + 5n 2 X 2 ) 



R 



5 3 



7X 



d 3 Re~ 2m y(KX,R) 



(25 + kA) 4 
A5 3 kX(25 + 5kX) 
(25 + nXf ' 



(37) 
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FIG. 4: The energy of the bbc ground state (in the heavier "baryonic" unit) as function of 
m re d/M re d. Three curves are generated by implementing three different approximation schemes. 
The dot-dashed line has the functional form E = — k — 2 be inferred from Q23[) . 

where C, £ are given in (J2EI), and 

y{\, R) = — [d 3 r e-^+r,) V R R-V Rri = ^ £(X, R) . (38) 

7T J O 

It is interesting to note that the contribution of the charm energy, which is previously 
computed in Born-Oppenheimer procedure, (|28jl . is also subsumed in (jHEJ) in a similar format. 
Besides this, (IHfiJ) also incorporates terms that have been neglected in Born-Oppenheimer 
approximation, such as y (see the comment in Footnote EJ) . 

The minimum of (j3*r]|) can be found by enforcing dE/d\\s = dE/dS\\ = 0. It is rather 
difficult to derive analytic expressions for these optima, hence we resort to numerical method 
to determine them. Subsequently, the energy of bbc ground state as function of k, juxtaposed 
with predictions made by two other approaches, is shown in Fig. 0] The optimized values 
of A and S as functions of k are shown in Fig. El 

As is expected, the energy predicted from this approach coincides with those from the 
other two in the k — > limit. The technical reason is easily traceable. Note all the integrals 
in (|3*rj)l simplify greatly in this limit, e.g., S, X m 1, and C, £ ~ A. Neglecting higher order 
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FIG. 5: Dependence of two optimized variational parameters A, 5 on the the mass ratio m ie &/M Te &. 



terms, Eq. (p3fi|) then reduces to 



E = -- + 8(8-1)+k 



-y + MA-2) 



0{k 2 



(39) 



The optima 5 = 1, A = 2 can be trivially inferred, and the corresponding energy is exactly 
the same as (}2*3*j) . which was first derived in the point-like diquark approximation. 

Fig. El illustrates some anticipated features of a bbc state. As k grows, this state starts 
to depart from the simple point-like diquark picture, and the effect of finite diquark size 
becomes increasingly important. It can be clearly observed that 5 ascends in a slower pace 
than A descends. This is compatible with the expectation that the impact of c on b is less 
important than the impact of b on c. 

One interesting observation from Fig.|3]is that Born-Oppenheimer approximation renders 
rather close prediction to that from the variational approach, virtually in all k range. The 
reason is perhaps that those terms dropped by Born-Oppenheimer procedure turn out to be 
insignificant in this case. In any rate, this approximation scheme is not expected to work so 
well when k gets large. It is worth mentioning that, when analyzing baryon mass spectrum 
from potential model approach, Fleck and Richard have also found this scheme yields rather 
accurate results j^j. They have attributed it to a lore that asserts Born-Oppenheimer 
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approximation works always better than expected. 

We end this section by pointing out an interesting finding. In the complicated expression 
for E, (|36|) . the last two terms nearly cancel with each other in virtually all the range of k, 
once the optimized values of A and 5 are used. We thus achieve a great simplification: 

A 2 X 2 

S-4-4- (40) 

This approximate formula works surprisingly well. It deviates from the actual one by 3% 
in maximum in the range K < 1. If one restricts to the smaller range k < 0.8, the error of 
this formula is less than 1 %. 

Without a deeper understanding, one might simply regard the success of (jjOj) as a fortu- 
itous coincidence. If taken seriously, it seems to indicate that the point-like diquark picture 
might be useful even at large k. It will yield the right answer, if one pretends that the color 
potential between two b quarks is —5/R, and the bb diquark perceived by c is equivalent to 
an antiquark carrying the color charge A. 

IV. PHENOMENOLOGY 

In this section, we will assemble the knowledge gleaned in the preceding section to esti- 
mate masses of various lowest-lying triply heavy baryons. We then compare our results with 
other work in literature, and discuss corresponding implications. 

It should be first realized that our predictions will be sensitive to the input of heavy 
quark masses. Therefore, it is important to specify an appropriate quark mass scheme to 
lessen arbitrariness. Since our working assumption is the weak-coupling regime, it is most 
consistent to express the heavy quark pole mass in terms of the masses of lowest-lying 
quarkonia, J/^f and T, assuming they are the weakly coupled system. At order a 2 s , we can 
write m c and m& as 



m b 



Mj,y 

2 

My 



1 

9 



1 + 



9 



(41) 



We will take the physical values Mj^ = 3.097 GeV, M T = 9.460 GeV as input. 

We are now at a position to express the masses of tripled-heavy baryons in perturbative 
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expansion. We start from flbcc- Using the result of (JUJ), we find 



mi, + 2m c + E 



2 



+ M J/9 + 



\M m + \M T 



8j 2(M J/9 + M r )_ 



+ Mj/* [1 + 0.273 a*(//)] . 



(42) 



In expressing the reduced mass, we simply replace m c with half of Mjm, and m& with half 
of Mt. This simplified procedure induces an error of order to the baryon mass, thus 
legitimate at present 0(a 2 s ) accuracy. 

The masses of baryons made of three identical quarks can be estimated in a similar 
manner. With the input from ()2U|). we infer the Q ccc mass to be 



Mr, 



3m c + E 



3M 



j/4 



0.925 
6 



3M 



j/i> 



1 + 0.068 a%(ji) 



(43) 



and the Qm mass can be obtained by making obvious replacement. 

For the bbc state, we have attempted three different approaches to estimate the binding 
energy. Since the one-step variational estimate is believed to be most reliable, we will adopt 
its prediction (though Born-Oppenheimer approximation yields a close result, as disclosed 
in Fig. |3J). First we need specify the value of k, the ratio of two reduced quark masses. It 
can be approximated as 

4M J/vI , 



K 



0.563 



Mj^ + 2 My 

and the error brought in by this procedure is assumed to be negligible. 

The optima can be determined numerically from (|3*Bj) by variational ansatz: 



(44) 



5 



1.137. 



A = 1.603. 



(45) 



with the corresponding energy 

E = -1.363 



-1.363 M red 



(2a t 



V 3 



(46) 



where we have inserted the Bohr energy of b quark in the last entity. 
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Piecing everything together, we obtain 



2m& + m c + E 
2 



+ M T + M j/ip + (2 - 1.363)M T 
+ M T [1 + 0.107 a^)]. 



a£(p) 



(47) 



So far we have treated each of triply-heavy baryons separately, it is not yet clear whether 
there is any connection among them. Interestingly, there is a mass convexity inequality 
relating different baryon states, which arises from general reasoning in QCD [21| . To our 
interest, such an inequality demands 



Ma 



(48) 



The underlying assumption of this theorem is universal interquark potential. Taking a s in 
(jUJ), an d (|4T|) to be equal, we readily verify that our predictions based on variational 
ansatz are indeed compatible with this QCD theorem. 
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There also exists another inequality, which relates the masses of baryons and mesons 
l^ . This is derived from the assumption that the quark-quark potential in a baryon is a half 
of the quark- ant iquark potential in a meson, which is de facto satisfied in Coulomb bound 
states. To our purpose, this inequality reads 



Mo 



hhc — 2 



+ M £ 



(49) 



To make a consistent examination of this relation, we need treat B c also as a weakly- 
coupled state, which is believed to be the case. Following preceding procedure, we can 
express the B c mass as 

4M J/q ,M r 



Mf 



M r | Mj/y | 



2 

M J/q , 



+ 



2 

M r 



M r + Mj,^ 



M J/q , + M T 



9 



[1 + 0.076 a s 2 (/x)] 



2 2 

One can promptly check that this inequality also holds in our case. 

A simple variant of (}4*9~|) is to specify all the quarks to be of a single flavor 



(50) 



221: 
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> 



3M 



j/i> 



(51) 



Our prediction in (|4^j) indeed respects this requirement. 
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TABLE I: Predictions for the masses of lowest-lying triply-heavy baryons from various work. All 
the masses are given in unit of GeV. In the entries for £lbcc and &bbc, the J p = ^ + and J = | + 
partners are not distinguished since the hyperfine splitting has been neglected. 





Bjorken [4] 


This work 


Vijande et al [24] 


^bcc 


8.200 ± 0.090 


7.98 ±0.07 




^ccc 


4.925 ± 0.090 


4.76 ± 0.06 


4.632 


Qbbb 


14.760 ±0.180 


14.37 ±0.08 




&bbc 


11.480 ±0.120 


11.19 ±0.08 





To make quantitative estimates for the baryon masses, we need specify at which scale the 
strong coupling constant should be evaluated. In principle, physical observables should be 
independent of the choice of /i, once the all-order perturbative expansion has been worked 
out. In practice, since what we have so far is only the leading order perturbative correction, 
our predictions are unavoidably sensitive to the choice of /i. To reduce the scale ambiguity 
optimally, we should take [i in proximity to the characteristic momentum transfer scale in 
a given QQQ state. 

It is an empirical fact that the typical momentum transfer scale in J/ip, B c and T is 
about 0.9, 1.2 and 1.5 GeV, respectively. One might expect that the corresponding scale in 
the QQQ states would be considerably lower than that in their quarkonium counterparts. 
Encouragingly, as we have learned in Sec. IIII| the effective color strength between a pair of 
quarks gets enhanced due to the presence of the third quark. As a result, the actual wave 
function is more compressed than naively expected. Therefore, it is not unreasonable to 
choose the scale for a QQQ state close to the one typically taken for its QQ counterpart. 
We assign \x = 1.2 GeV in the mass formula for flbcc, &bbb and flbbc, with a corresponding 
a s = 0.43; for fl ccc , we take /i = 0.9 GeV, with a s = 0.59. To compensate for our ignorance 
in uncalculated higher order corrections, we estimate the uncertainty in each mass prediction 
to be the leading 0(a1) correction multiplied by another factor of a s . 

Our predictions to the masses of various QQQ ground states, together with those made 
by other authors 0, |3| , which employ some phenomenological confinement potentials, are 
compiled in Table HI The apparent discrepancy between the predictions of the Q ccc mass 
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by Bjorken and by Vijande et al, which is as large as 300 MeV, might reflect the large 
uncertainty inherent in phenomenological approaches 4 . In contrast, our predictions are 
based on the perturbation theory, being systematically improvable, suffer less arbitrariness. 

It can be readily recognized that Bjorken's predictions are systematically higher than 
ours. Note the variational method by default underestimates the binding energy, and a 
more accurate weak-coupling analysis will predict even lower masses for QQQ ground states, 
hence further enlarging this disagreement. 

Note that the 0{af) corrections in (JUJ), (SSI) and (|37)l are all positive, so as we lower 
down fi, which is meant to be the characteristic momentum scale, our predictions will shift 
upwards, getting close to Bjorken's predictions. When \i descends further and becomes 
comparable with Aqcd> our method breaks down and one enters the strong-coupling regime. 
In a sense, Bjorken's results can be considered as arising from a strong-coupling analysis. 

The future experiments and lattice QCD simulations will decide which prediction is closer 
to the reality, consequently nature of the QQQ ground states may be disclosed. 

V. SUMMARY AND OUTLOOK 

The theme of this work is to estimate the masses of various lowest-lying triply heavy 
baryon states, with the assumption that they are weakly-coupled system, analogous to T, 
B c and J/ip. To achieve this, it is crucial to make a sound estimate for the binding energy 
of a nonrelativistic three heavy quark system, which is bound by short-distance interquark 
potentials that are organized by powers of a s and 1/M. Due to our incapability of rigorously 
solving 3-body problem, we have invoked the variational method as an approximation scheme 
to analyze various QQQ ground states. As a first step, we have estimated the most important 
piece, i.e., the 0(a 2 s ) contribution to the binding energy, with only the tree-level static 
potential incorporated. 

For the variational method to be accurate, it is important to choose a reasonable form of 
trial states. In view of this, different triply-heavy baryon states, the bcc, ccc (bbb) and bbc 
states, have been analyzed separately, each supplied with a different trial wave function mo- 
tivated by the symmetry consideration and the presence of hierarchy m& 3> m c . Inspired by 

4 Note the very low fl ccc mass predicted by Vijande et al violates the mass inequality (|51|l . 
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the similarity between our baryonic system and the three-body atomic system, some guid- 
ances have been taken from the familiar textbook treatment of helium atom and the ionized 
hydrogen molecule. Among various QQQ states, the bbc state is the most interesting one but 
most challenging to analyze. We have carried out a detailed study on this state, employing 
several different approaches. The implications of different approaches are elucidated, and in 
particular the relevance of the compact diquark picture has been discussed. 

Masses of various QQQ ground states derived from our formalism are compatible with 
those well-known mass inequalities in QCD. Our quantitative predictions, which is based 
on a weak-coupling treatment, are systematically lower than Bjorken's, which may instead 
be viewed as resulting from a strong-coupling analysis. It leaves for future experiments 
and lattice QCD simulation to decide the nature of the lowest-lying triply heavy baryons, 
whether to be weakly coupled or strongly coupled. 

Besides the ability to estimate the masses, variational analysis also allows us to have a 
reasonable knowledge about the quark wave functions. This will be useful, for example, in 
estimating the hyperfme splittings in the bcc and bbc states. One important byproduct of 
this knowledge is that a reasonable value of the wave function at the origin can be inferred. 
Like in a heavy quarkonium, this quantity is one of the basic characteristics of a triply-heavy 
baryon state, and is of phenomenological interest. For instance, this value is a crucial input 
for reliably estimating the fragmentation function for the QQQ states f|. 

It will be useful if alternative methods that have been developed to analyze 3-body 
problem, e.g., the formalism of hyperspherical expansion [3], are employed to check the 
accuracy of our results. 

One apparent improvement on this work is in prospect. As stressed several times before, 
we have only incorporated the tree-level static potential in this work, so that our predictions 
suffer from considerable scale dependence. The perturbative matching calculation for the 
one-loop static potential and tree-level spin-dependent potentials is straightforward. It will 
be useful to implement their contributions into our variational framework. 
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